Computer method and apparatus for optimized controller in a non-linear process

ABSTRACT

A first principles, steady state model of a desired polymer process is applied with a non-linear optimizer to a linear controller. Model process gains and optimal target values for controller variables result. These results are utilized by a multivariable linear controller to achieve nonlinear control of the subject process. Preferably the nonlinear optimizer is DMO/SQP. The steady state model is produced by Polymers Plus and the linear controller is DMCplus, all of Aspen Technology, Inc. in Cambridge Mass.

RELATED APPLICATION

[0001] This application claims the benefit of Provisional PatentApplication No. 60/171,799, filed Dec. 22, 1999, the entire teachings ofwhich are incorporated herein by reference.

BACKGROUND OF THE INVENTION

[0002] The polymer process is a complex nonlinear process. There are,therefore, many types of processes developed by different manufacturers.The differences within a single product type, such as polyethylene,include process configuration (e.g. tubular reactors, stirred tankreactors, loop reactors), reaction medium (e.g. gas phase, slurry,solution), catalyst types (Ziegler-Natta, peroxide, chromium, vanadium,and metallocene), reaction pressure and reaction temperature. As aconsequence, these polymer processes exhibit significantly differentnonlinear effects upon product properties.

[0003] For most polymer processes, the operating characteristics involvemaking one type of product for a period of time to satisfy a productorder and then changing operating conditions to make another producttype for a new demand. Typically, product types are characterized bybulk properties such as Melt Index and Density, which indicate how theproduct will behave when it is moulded or blown into a film. There aremany other variations of these measurements, as well as other visual andperformance properties, such as color and fish eyes, that are much moredifficult to predict and control. These differences in design andcharacterization vary even more across products such as polypropylene,polystyrene, polycarbonates, nylon, etc.

[0004] Historically, it has been a challenge to control industrialpolymer processes. Currently, the standard practice is to use neuralnetwork regression to identify process gains needed to adapt amultivariable linear controller in order to achieve a kind of nonlinearcontrol. Aspen IQ™ and DMCplus™ (both by Aspen Technology, Inc. ofCambridge, Mass.) are examples of such a neural network program andlinear process controller, respectively. The DMCplus linear models arebased on linearized models around the nominal operating point. Thecurrent model gains are used by DMCplus for calculation of the gainmultipliers. However, this approach has proven to be time consuming,manpower intensive and costly.

SUMMARY OF THE INVENTION

[0005] The present invention provides a solution to the foregoingproblems in process control in the prior art. In particular, the presentinvention provides a computer method and apparatus which enables amultivariable, process controller to achieve non-linear control. In apreferred embodiment, the present invention utilizes the rigorous,non-linear model of the process at steady state as generated by PolymerPlus® (a software product by Aspen Technology, Inc. of Cambridge, Mass.)to optimize the controller.

[0006] Hence, in accordance with one aspect of the present invention, anonlinear optimizer solves a first principles, steady state processmodel and calculates process gains and optimal targets for themultivariable controller. The first principles, rigorous, mechanisticPolymers Plus models handle the issue of process non-linearity derivedfrom kinetics, thermodynamics and process configuration. These modelsare valid across a wide operating range, extrapolate well, capture theprocess non-linearity and require only minimal amounts of process data.Based on this approach, the current process gains for eachIndependent/Dependent model can be easily obtained from the partialderivatives of the corresponding first principles Polymers Plus model.

[0007] In the preferred embodiment, the optimizer calculates the optimaltargets for the Manipulated Variables (MVs) and Controlled Variables(CVs) of the DMCplus controller, replacing the internal Linear Program(LP) optimizer that is, based on the current process gains. This way,the DMCplus controller follows a consistent set of targets and does notchange its direction due to process gain changes. It is noted that theDMCplus controller still uses the current model gains (based on thecurrent gain multipliers) to calculate the control-move plan so thatcontroller stability is preserved.

[0008] To that end, computer apparatus embodying the present inventioncomprises (a) a controller for determining and adjusting manipulatedvariable values for controlling a subject non-linear manufacturingprocess, and (b) an optimizer coupled to the controller for updating thelinear model of the controller. The controller employs a dynamic linearmodel for modeling the effect that would result in the subjectmanufacturing process with a step change in manipulated variable values.As the subject non-linear manufacturing process transitions from oneoperating point to another, in a high degree of non-linearity betweenmanipulated variables and controlled variables of the subject process,the optimizer updates the linear models of the controller. The optimizerutilizes a non-linear model of the subject process for determiningtarget values of the controlled variables. The controlled variables areindicative of physical properties of the subject process.

[0009] In accordance with another aspect of the present invention, thereis a source of sensor measured variables for representing the measurablephysical properties and hence controlled variables of the subjectprocess. The non-linear model of the optimizer determines gains betweenthe manipulated variables and the sensor measured controlled variables.As such, the optimizer gain adapts the linear model of the controllerwith the determined gains.

[0010] In accordance with another aspect of the invention, thenon-linear model of the optimizer is a rigorous, first principles,non-linear model. Further, the optimizer and its non-linear model isexecuted as frequently as the controller.

[0011] The present invention method for controlling a non-linearmanufacturing (e.g., polymer) process thus includes thecomputer-implemented steps of:

[0012] (i) utilizing a linear model, modeling effect that would resultin a subject manufacturing process with a step change in manipulatedvariable values used for controlling said process;

[0013] (ii) using a non-linear model of the subject process, determiningtarget values of the controlled variables indicative of physicalproperties of the subject process; and

[0014] (iii) updating the linear model as the subject processtransitions; from one operating point to another, in a high degree ofnon-linearity between the manipulated variables and controlled variablesof the subject process.

[0015] In particular, the invention method uses the non-linear model ofthe subject process to update the process gains (between the manipulatedvariables and the controlled variables) for the linear model.

BRIEF DESCRIPTION OF THE DRAWINGS

[0016] The foregoing and other objects, features and advantages of theinvention will be apparent from the following more particulardescription of preferred embodiments of the invention, as illustrated inthe accompanying drawings in which like reference characters refer tothe same parts throughout the different views. The drawings are notnecessarily to scale, emphasis instead being placed upon illustratingthe principles of the invention.

[0017]FIG. 1 is a schematic drawing of a manufacturing process with anoptimized controller of the present invention.

[0018]FIG. 2 is a block diagram of the controller of FIG. 1 asimplemented in a computer system.

[0019]FIG. 3 is a block schematic of a reactor model employed in theoptimizer to the controller of FIG. 2.

DETAILED DESCRIPTION OF THE INVENTION

[0020] Illustrated in FIG. 1 is a manufacturing plant for carrying out anon-linear process such as a chemical or polymer process. The processingplant is formed of a series of vessels (e.g., holding tanks, feed tanks,catalyst feeds), reactors (including mixing tanks, etc.) and pumps (orcompressors) connected by various conduits and pipes. Sensors 17 sensestemperature, volume, pressure, compositions and other physicalproperties at various points in the process and provide sensor data to asoftware or computer-implemented controller 19. The controller 19utilizes the sensor data to maintain setting for controlled variablessuch as temperature, pressure, compositions and product properties byadjusting manipulated variables such as feed rates, flowrate,temperature of the vessels, reactor 13 and pumps/compressors 15.Controller 19 physically effects adjustment of the manipulated variablesthrough actuators 21 coupled to respective valves 23 of the vessels 11,reactors 13 and pumps/compressors 15.

[0021] In particular, the controller 19 employs a linear dynamic modelof the manufacturing process. These linear models relate the dynamicresponses of controlled variables to manipulated variables in terms ofprocess gains, response time, and dead time. The sensor data definevalues of control variables in the model equations. The model predictshow controlled variables will change with respect to step changes inmanipulated variables.

[0022] An optimizer 25 uses the sensor data and an internal model of themanufacturing process to provide target values for the controlledvariables. In the preferred embodiment, commonly known first principlesequations for thermodynamics, kinetics, and heat and mass balancesdefine the model. As a result, the model is capable of predicting thenonlinear relationship (or gains) between the controlled variables andthe manipulated variables based on first principles.

[0023] In the case of a highly non-linear process being controlled, thelinear model of the controller 19 is problematic. A common practice inthe prior art is to update the process gains of the linear model in apre-programmed manner based on process experiences. Applicants havediscovered an improved non-linear optimizer 25 for such a controller 19.Specifically, the optimizer 25 of the present invention employs anon-linear first principles model which is suitable for frequentlycalculating target controlled variable values and non-linear gainsbetween the controlled variables and manipulated variables. Frequency ofoptimizer 25 calculations and operation is about every 2-3 minutes orabout the same rate as the controller 19 determinations are made. Thenon-linear gains computed by the optimizer 25 are used continuously toupdate the process gains of the linear model of the controller 19. Incontrast, typical optimizers used with process control of the prior artperform calculations only for target control variables at about onceevery four hours or so and hence not at a rate helpful to the controller19 of a highly non-linear process.

[0024] Referring now to FIG. 2 is an illustration of the preferredembodiment in a computer system 31. Generally, the computer system 31includes a digital processor 33 which hosts and executes the inventioncontroller 19 and optimizer 25 in working memory. Input 35 to digitalprocessor 33 is from sensors 17, another software program, anothercomputer, input devices (e.g., keyboard, mouse, etc.) and the like.Output 37 (i.e., the controlled variable settings determined bycontroller 19 and gain adapted by optimizer 25) is provided to actuators21, another computer, storage memory, another software program and/oroutput devices (e.g., display monitor, etc.) and the like. It isunderstood that computer system 31 may be linked by appropriatecommunication links to a local area network, external network (i.e., theInternet) or similar such networks for sharing or distributing input andoutput data.

[0025] In FIG. 2, the controller 19 is preferably DMC Plus by AspenTechnology, Inc. of Cambridge, Mass. Optimizer 25 is preferablyDMO/SQP®) also of Aspen Technology, Inc. Other similar non-linearoptimizers are suitable. In the preferred embodiment, Polymers Plus(module 27 in FIG. 2) of Aspen Technology Inc. supplies the internalrigorous, non-linear model 29 to optimizer 25. That is, optimizer 25uses a steady state process model 29 based on first principles,rigorous, mechanistic Polymers Plus models.

[0026] By way of the below example, the preferred embodiment isillustrated around a fluidized bed, gas phase polyethylene process. Thenon-linear optimizer 25 solves the non-linear steady state model 29 fromPolymer Plus module 27. To accomplish this, the optimizer 25 employs a:

[0027] Sparsity file,

[0028] Nonlinear steady state model 29, and

[0029] Objective function

[0030] The non-linear model 29 is, preferably, formulated inopen-equation form. In this case, optimizer 25 supplies values to themodel 29 for all the variables of interest, including manipulatedvariables and controlled variables. The model 29 returns to optimizer 25the values of the constraint equation residuals, and as many of theJacobian elements (partial derivatives) as it can calculate.

[0031] If an open-equation model is not available, a closed-form modelmay also be used. In that case, optimizer 25 supplies values to themodel 29 for its ‘input’ variables, and the model 29 returns values ofits ‘output’ variables.

[0032] The nonlinear model 29 user interface package allows users to:

[0033] Set a scenario—initial values of controlled variables andmanipulated variables

[0034] Set constraints—targets or upper and lower limits on controlledvariables and manipulated variables

[0035] Set the objective function—costs on controlled variables andmanipulated variables

[0036] Run a simulation of the scenario

[0037] View the calculated trajectory of the controlled variable targetsand manipulated variable targets over the nonlinear simulation interval.

[0038] As such, Optimizer 25 (via the nonlinear model 29) supplies thefollowing to the DMCplus controller 19:

[0039] Controlled variable targets or upper and lower limits

[0040] Manipulated variable targets

[0041] Model gains (i.e. derivatives of controlled variables withrespect to manipulated variables, derivatives of controlled variableswith respect to Feed Forward variables) for the current operating point

[0042] Other model variables (such as residence time)

[0043] Example Reactor Model

[0044] The nonlinear model 29 incorporates all of the equipment ofinterest required to represent the polyethylene process—the reactor,heat exchanger, compressor, mixers, splitters, component separator,valves, etc. Details of the reactor model 49 are given below forillustrative purposes.

[0045] As shown in FIG. 3, the reactor 49 is modeled as a steady-stateCSTR (continuous stir tank reactor), with a vapor phase 45 and a liquidphase 47 in equilibrium. Two feeds—a gas feed 41 and a catalyst feed43—are required. The two products are (1) a vapor stream 51, and (2) aliquid stream 53 containing the polymer.

[0046] The reactor model 49 is composed of the following nonlinearequations derived from first principles of the Polymers Plus module 27:

[0047] Component material balances

[0048] Total material balance

[0049] Pressure balances

[0050] Energy balance

[0051] Vapor enthalpy calculation

[0052] Liquid enthalpy calculation

[0053] Vapor-liquid equilibrium

[0054] Reactor volumes

[0055] Component reaction rates

[0056] Polymer attribute calculations

[0057] Catalyst attribute calculations

[0058] The equations require calculations on the stream enthalpy, streamdensity, reaction kinetics, vapor-liquid equilibrium K-values, and thepolymerization reaction kinetics. Such calculations are performed bysubroutines in Polymers Plus module 27.

[0059] Component Slate

[0060] The component slate for the reactor model 49 consists of CatCatalyst Cocat Cocatalyst C2 Ethylene C4 Butene H2 Hydrogen HDPE Polymer

[0061] Variables

[0062] Streams into and out of the reactor model 49 have a standardformat and units of measurement to enable automatic connection to otherequipment models. Conversion to internal model units of measurementtakes place inside the model 49. The stream variables are: Gas feed(41): F_(gas) Gas feed flow Klbmol/hr Z_(l) Gas feed component molefractions mole fraction T_(gas) Gas feed temperature deg F P_(gas) Gasfeed pressure PSIG H_(gas) Gas feed enthalpy KBTU/lbmol Catalyst Feed(43): F_(cat) Catalyst flow Klbmol/hr W_(l) Catalyst component molefractions mole fraction T_(cat) Catalyst temperature deg F P_(cat)Catalyst pressure PSIG H_(cat) Catalyst enthalpy KBTU/lbmol VaporProduct (51): F_(V) Vapor flow Klbmol/hr y_(i) Vapor component molefractions mole fraction T_(V) Vapor temperature deg F P_(V) Vaporpressure PSIG H_(V) Vapor enthalpy KBTU/lbmol Liquid Product (53): F_(L)Liquid flow Klbmol/hr x_(i) Liquid component mole fractions molefraction T_(L) Liquid temperature deg F P_(L) Liquid pressure PSIG H_(L)Liquid enthalpy KBTU/lbmol Reactor Model 49 variables: Q Heat added toreactor MBTU/hr Level Liquid level in reactor meters V_(L) Volume ofliquid cubic meters V_(V) Volume of vapor cubic meters Rho_(L) Liquiddensity kgmol/cum Rho_(V) Vapor density kgmol/cum ResTime_(L) Liquidresidence time hours ResTime_(V) Vapor residence time hours Pol_(out)Polymer flow out kg/sec R_(l) Component reaction rate (−ve =consumption) kgmol/cum-sec R_(SZMoM) Reaction rate for zeroth moment ofbulk polymer kgmol/cum-sec R_(SSFLOWseg) Reaction rate for first momentof bulk polymer kgmol/cum-sec R_(SSMOM) Reaction rate for second momentof bulk polymer kgmol/cum-sec R_(LSEFLOWseg) Reaction rate for zerothmoment of live polymer kgmol/cum-sec R_(LSSFLOWseg) Reaction rate forfirst moment of live polymer kgmol/cum-sec SZMOM Zeroth moment of bulkpolymer gmol/kg polymer SSFLOW_(seg) First moment of bulk polymer, persegment gmol/kg polymer SSMOM Second moment of bulk polymer kgmol/kgpolymer LSEFLOW_(seg) Zeroth moment of live polymer, per segmentmilli-gmol/kg polymer LSSFLOW_(seg) First moment of live polymer, persegment milli-gmol/kg polymer R_(CPS) Reaction rate for catalystpotential sites kgmol/cum-sec R_(CVS) Reaction rate for catalyst vacantsites kgmol/cum-sec R_(CDS) Reaction rate for catalyst dead siteskgmol/cum-sec Cat_(out) Catalyst flow out kg/sec CPS_(in) Catalystpotential site concentration-catalyst feed milli-gmol/kg catalystCVS_(in) Catalyst vacant site concentration-catalyst feed milli-gmol/kgcatalyst CDS_(in) Catalyst dead site concentration-catalyst feedmilli-gmol/kg catalyst CPSFLOW Catalyst potential siteconcentration-liquid product milli-gmol/kg catalyst CVSFLOW Catalystvacant site concentration-liquid product milli-gmol/kg catalyst CDSFLOWCatalyst dead site concentration-liquid product milli-gmol/kg catalystMWW Weight-average degree of polymerization MWN Number average degree ofpolymerization MI Polymer melt index MIBias Melt index offset frommeasured Frac Fraction comonomer in polymer mole fraction Dens Polymerdensity Dbias Polymer density offset from measured Internal variablesK_(l) Component K-value MW_(l) Component molecular weight g/gmolParameters The following variables have fixed values V Reactor volumecubic meters Area Reactor cross-section area, liquid section squaremeters A, B, C, D Constants in melt index equation E, F, G Constants inpolymer density equation E3 1000.0 E6 1.0E+06 Inputs to subroutineZNMECH: NSITES Total number of site types 1 NCAT Number of catalysts 1NCCAT Number of cocatalysts 1 NMOM Number of monomers 2 NSEG Number ofsegments 2 NPOL Number of polymers 1 AKO(nrx) Pre-exponential factorsEACT(nrx) Activation energies ORD(nrx) Reaction order TREF(nrx)Reference temperature 1.0E+35 Conc(ncpt) component concentrations =x_(i).Rho_(L) kgmol/cum CPS catalyst site concentration =CPSFLOW.Rho_(L).x_(cat).MW_(cat)/E6 kgmol/cum CVS catalyst siteconcentration = CVSFLOW.Rho_(L).x_(cat).MW_(cat)/E6 kgmol.cum Mu0(seg)live moment concentration = LSEFLOW_(seg).Rho_(L).x_(pol).MW_(pol)/E6kgmol/cum Mu1(seg) live moment concentration =LSSFLOW_(seg).Rho_(L).x_(pol).MW_(pol)/E6 kgmol/cum Lam0 dead momentconcentration = SZMOM.Rho_(L).x_(pol).MW_(pol)/E3 kgmol/cum Lam1(seg)dead moment concentration = SSFLOW_(seg).Rho_(L).x_(pol).MW_(pol)/E3kgmol/cum Lam2 dead moment concentration =SSMOM.Rho_(L).x_(pol).MW_(pol) kgmol/cum Conversion factors: Kg2LbConvert kilograms to pounds   2.2046 Sec2Hr Convert seconds to hours3600.0 Subroutines (of Polymers Plus module 27) DENSITY Stream densityENTHLP Stream enthalpy KVALUE Vapor-liquid equilibrium Kvalues ZNMECHComponent reaction rates, catalyst attributes, and polymer attributesEquations (solved by the non-linear model 29) Component materialbalances: Catalyst: F_(cat).w_(i) − F_(L).x_(i) +Sec2Hr.Kg2Lb.R_(i).V_(L)/E3 = 0 Klbmol/hr (1) Cocatalyst: F_(cat).w_(i)− F_(L).x_(i) + Sec2Hr.Kg2Lb.R_(i).V_(L)/E3 = 0 Klbmol/hr (2) Ethylene:F_(gas).z_(l) − F_(v).y_(l) −F_(L).x_(i) + Sec2Hr.Kg2Lb.R_(l).V_(L)E3 =0 Klbmol/hr (3) Butene: F_(gas).z_(l) − F_(v).y_(l) −F_(L).x_(i) +Sec2Hr.Kg2Lb.R_(l).V_(L)E3 = 0 Klbmol/hr (4) Hydrogen: F_(gas).z_(l) −F_(v).y_(l) −F_(L).x_(i) + Sec2Hr.Kg2Lb.R_(l).V_(L)E3 = 0 Klbmol/hr (5)Polymer: − F_(L)x_(l) + Sec2Hr.Kg2Lb.R_(l).V_(L).E3 = 0 Klbmol/hr (6)Total material balance: F_(gas) + F_(cat) − F_(L) +Sec2Hr.Kg2Lb.ΣR_(l).V_(L)E3 = 0 Klbmol/hr (7) Temperature balance: T_(V)− T_(L) = 0 deg F (8) Pressure balances: P_(cat) − P_(gas) = 0 PSIG (9)P_(V) − P_(gas) = 0 PSIG (10) P_(L) − P_(gas) = 0 PSIG (11) Energybalance: F_(gas).H_(gas) + F_(cat).H_(cat) − F_(V).H_(V) − F_(L).H_(L)+Q = 0 MBTU/hr (12) Vapor enthalpy calculation: H_(V) − Enthlp(‘V’,y_(i), T_(V), P_(V)) = 0 KBTU/lbmol (13) Liquid enthalpy calculation:H_(L) − Enthlp(‘L’, x_(l), T_(L), P_(L)) = 0 KBTU/lbmol (14)Vapor-liquid equilibrium: Ethylene, Butene, Hydrogen: y_(l) −K_(l).x_(l) = 0 mole fraction (15-17) Flash equation Σy_(l) − Σx_(l) = 0(18) Kvalue calculation K_(l) − K_(value)(x_(i), y_(l), T_(L), P_(L)) =0 Reactor volumes: Liquid volume: V_(L) − Level.Area = 0 cubic meters(19) Liquid density: Rho_(L) − Density(‘L’, x_(l), T_(L), P_(L)) = 0kgmol/cum (20) Liquid residence time: ResTime_(L) −V_(L).Rho_(L).Kg2Lb/F_(L) = 0 hours (21) Vapor volume: V − V_(L) − V_(V)= 0 cubic meters (22) Vapor density: Rho_(V) − Density(‘V’, y_(l),T_(V), P_(V)) = 0 kgmol/cum (23) Vapor residence time: ResTime_(V) −V_(V).Rho_(V).Kg2Lb/F_(V) = 0 hours (24) Component reaction rates: R_(l)− ZNMECH(T_(L), P_(L), x_(l), MW_(l), Mu_(l), Lam_(l), CPS, CVS, CIS) =0 kgmol/cum-sec (25-30) Polymer attributes: Attribute rates: R_(SZMOM) −ZNMECH(T_(L), P_(L), x_(l), MW_(i), Mu_(l), Lam_(i), CPS, CVS, CIS) = 0kgmol/cum-sec (31) R_(SSFLOWseg) − ZNMECH(T_(L), P_(L), x_(i), MW_(l),Mu_(l), Lam_(i), CPS, CVS, CIS) = 0 kgmol/cum-sec (32-33) R_(SSMOM) −ZNMECH(T_(L), P_(L), x_(l), MW_(i), Mu_(l), Lam_(l), CPS, CVS, CIS) = 0kgmol/cum-sec (34) R_(LSEFLOWseg) − ZNMECH(T_(L), P_(L), x_(l), MW_(i),Mu_(l), Lam_(l), CPS, CVS, CIS) = 0 kgmol/cum-sec (35-36) R_(LSSFLOWseg)− ZNMECH(T_(L), P_(L), x_(l), MW_(l), Mu_(l), Lam_(l), CPS, CVS, CIS) =0 kgmol/cum-sec (37-38) Polymer flow out Pol_(out) −F_(L).x_(pol).MW_(pol).E3/Sec2Hr.Kg2Lb = 0 kg/sec (39) Zeroth moment ofbulk polymer, per site (1 site) SZMOM.Pol_(out)/E3 − R_(SZMOM).V_(L) = 0kgmol/sec (40) First moment of bulk polymer, per segment (2 segments),per site (1 site) SSFLOW_(seg).Pol_(out)/E3 − R_(SSFLOWseg).V_(L) = 0kgmol/sec (41, 42) Second moment of bulk polymer, per site (1 site)SSMOM.Pol_(out)/E3 − R_(SSMOM).V_(L) = 0 kgmol/sec (43) Zeroth moment oflive polymer, per segment (2 segments), per site (1 site)LSEFLOW_(seg).Pol_(out)/E6 − R_(LSEFLOWseg).V_(L) = 0 kgmol/sec (44, 45)First moment of live polymer, per segment (2 segments), per site (1site) LSSFLOW_(seg).Pol_(out)/E6 − R_(LSSFLOWseg).V_(L) = 0 kgmol/sec(46, 47) Melt index MWW − SSMOM/ΣSSFLOW_(seg) = 0 weight-average (48)MWN − ΣSSFLOW_(seg)/SZMOM = 0 number-average (49) MI − A.MWW^(B) −C.MWN^(D) + MIBias = 0 (50) Density Frac − SSFLOW_(C4)/ΣSSFLOW_(seg) = 0mol fraction (51) Dens − E + F.Frac^(G) + Dbias = 0 (52) Catalystattributes: Attribute rates R_(CPS) − ZNMBCH(T_(L), P_(L), x_(l),MW_(l), Mu_(l), Lam_(i), CPS, CVS, CIS) = 0 kgmol/cum-sec (53) R_(CVS) −ZNMECH(T_(L), P_(L), x_(l), MW_(l), Mu_(i), Lam_(i), CPS, CVS, CIS) = 0kgmol/cum-sec (54) R_(CDS) − ZNMECH(T_(L), P_(L), x_(i), MW_(i), Mu_(l),Lam_(i), CPS, CVS, CIS) = 0 kgmol/cum-sec (55) Catalyst flow outCat_(out) − F_(L).x_(cat).MW_(cat).E3/Sec2Hr.Kg2Lb = 0 kg/sec (56)Potential sites, per catalyst (1 catalyst) CPSFLOW.Cat_(out)/E6 −CPS_(in).F_(cat).w_(cat).MW_(cat)/Sec2Hr.Kg2Lb.E3 − R_(CPS).V_(L) = 0kgmol/sec (57) Vacant sites, per site (1 site) CVSFLOW.Cat_(out)/E6 −CVS_(in).F_(cat).w_(cat).MW_(cat)/Sec2Hr.Kg2Lb.E3 − R_(CVS).V_(L) = 0kgmol/sec (58) Dead sites CDSFLOW.Cat_(out)/E6 −CDS_(in).F_(cat).w_(cat).MW_(cat)/Sec2Hr.Kg2Lb.E3 − R_(CDS).V_(L) = 0kgmol/sec (59)

[0063] Analysis

[0064] Fixed variables are those having values specified to define theproblem. The number of equations must equal the number of calculatedvariables.

Number of Variables

[0065] Description Number Fixed Calculated Gas feed - 3 components, F,T, P, H  7  6  1 (P) Catalyst feed - 2 components, F, T, P, H  6  5  1(P) Vapor product - 3 components, F, T, P, H  7  1 (P)  6 Liquidproduct - 6 components, F, T, P, H 10  1 (T)  9 Heat added, Q  1  0  1Liquid Volume, Level, Residence time  7  1  6 Component reaction rates 6  0  6 Polymer attributes, rates 17  0 17 Melt index  4  1  3 Density 3  1  2 Catalyst attributes, rates 10 3  7 Total 78 19 59

Number of Variables

[0066] Description Number Component material balances  6 Total materialbalance  1 Temperature balance  1 Pressure balances  3 Energy balance,product enthalpies  3 Vapor-liquid equilibrium, 3 components + flash  4Reactor volume, Residence time  6 Component reaction rates  6 Polymerattributes, rates 17 Melt index  3 Density  2 Catalyst attributes, rates7 Total 59

[0067] In the preferred embodiment, the reactor model 49 is formulatedto fit into pre-defined stream conventions:

[0068] Mole fractions must sum to 1.0

[0069] A stream contains mole fractions, flow, temperature, pressure,and enthalpy

[0070] In Polymers Plus (module 27) formulation, the following variablesare included in the streams:

[0071] Catalyst feed stream 43 CPS_(in) Catalyst potential siteconcentration - catalyst feed milli-gmol/kg catalyst CVS_(in) Catalystvacant site concentration - catalyst feed milli-gmol/kg catalystCDS_(in) Catalyst dead site concentration - catalyst feed milli-gmol/kgcatalyst Liquid product stream 53 SZMOM Zeroth moment of bulk polymergmol/kg polymer SSFLOW_(seg) First moment of bulk polymer, per segmentgmol/kg polymer SSMOM Second moment of bulk polymer kgmol/kg polymerLSEFLOW_(seg) Zeroth moment of live polymer, per segment milli-gmol/kgpolymer LSSFLOW_(seg) First moment of live polymer, per segmentmilli-gmol/kg polymer CPSFLOW Catalyst potential site concentration -liquid product milli-gmol/kg catalyst CVSFLOW Catalyst vacant siteconcentration - liquid product milli-gmol/kg catalyst CDSFLOW Catalystdead site concentration - liquid product milli-gmol/kg catalyst

[0072] The manipulated variables of the reactor 49 include the catalystfeed 43 variables and gas feed 44 variables as solved by subroutineZNMECH. Thus, optimizer 25 forwards, in pertinent part, the output ofsubroutine ZNMECH to DMCPlus controller 19. In response, DMCPluscontroller 19 uses the subroutine outputs to update the process gains inits dynamic linear model of reactor 49.

[0073] The corresponding dynamic linear models of reactor 49 employed bythe DMC controller 19 are of the form

βΔu(k)=e(k+1)

[0074] where β is the reactor's dynamic matrix having m columns of thereactor's step function appropriately shifted down in order. Δu(k) is anm-dimensional vector of control moves. e(k+1) is the projected errorvector. See B. A. Ogunnaike and W. H. Ray, “Process dynamics, modeling,and control,” Chapter 27, pp. 1000-100.7, Oxford University Press 1994.

[0075] According to the foregoing example and description, the presentinvention utilizes first principles kinetics, thermodynamics and heatand mass balances to develop the non-linear relationships between themanipulated variables and controlled variables. Based on thesenon-linear relationships (non-linear model 29), optimizer 25 effectivelyand timely gain adapts the dynamic linear model for use in controller 19for the subject manufacturing process. As such, the preferred embodimentprovides a novel combination of (i) Polymers Plus rigorous, non-linearmodels 19, (ii) a non-linear optimizer 25 and (iii) DMCPlus (or similardynamic linear model based) controller 19 as heretofore unapplied toprocess control.

[0076] While this invention has been particularly shown and describedwith references to preferred embodiments thereof, it will be understoodby those skilled in the art that various changes in form and details maybe made therein without departing from the scope of the inventionencompassed by the appended claims.

What is claimed is:
 1. Computer apparatus for controlling a non-linearmanufacturing process comprising: a controller for determining andadjusting manipulated variable values for controlling a subjectnon-linear manufacturing process, the controller employing a linearmodel for modeling effect that would result in the subject manufacturingprocess with a step change in manipulated variable values; and anoptimizer coupled to the controller for updating the linear model of thecontroller, as the subject non-linear manufacturing process transitionsfrom one operating point to another, in a high degree of non-linearitybetween manipulated variables and controlled variables of the subjectprocess, the controlled variables being indicative of physicalproperties of the subject process, the optimizer utilizing a non-linearmodel of the subject process for updating the linear model anddetermining target values of the controlled variables.
 2. Computerapparatus as claimed in claim 1 wherein the controller employs dynamiclinear models.
 3. Computer apparatus as claimed in claim 1 furthercomprising a source of sensor measured variables for representing thecontrolled variables of the subject process.
 4. Computer apparatus asclaimed in claim 3 wherein the non-linear model of the optimizerdetermines gains between the manipulated variables and the sensormeasured controlled variables, and the optimizer gain adapts the linearmodel of the controller with the determined gains.
 5. Computer apparatusas claimed in claim 1 wherein the optimizer is executed as frequently asthe controller.
 6. Computer apparatus as claimed in claim 1 wherein thenon-linear model of the optimizer is a rigorous, first principles,non-linear model.
 7. Computer apparatus as claimed in claim 1 whereinthe subject non-linear manufacturing process is a polymer process.
 8. Amethod for controlling a non-linear manufacturing process comprising thecomputer implemented steps of: utilizing a linear model, modeling effectthat would result in a subject manufacturing process with a step changein manipulated variable values used for controlling said process; usinga non-linear model of the subject process, (i) updating the linear modelas the subject process transitions from one operating point to another,in a high degree of non-linearity between the manipulated variables andcontrolled variables indicative of physical properties of the subjectprocess, and (ii) determining target values of the controlled variables.9. A method as claimed in claim 8 wherein the step of utilizing a linearmodel includes dynamic linear models.
 10. A method as claimed in claim 8further comprising the step of providing sensor measured variables forrepresenting various controlled variables of the subject process.
 11. Amethod as claimed in claim 10 further comprising the step of determiningprocess gains between the manipulated variables and the sensor-measuredcontrolled variables using the non-linear model; and the step ofupdating the linear model includes gain adapting the linear model withthe determined gains.
 12. A method as claimed in claim 8 wherein thestep of updating is executed as frequently as the subject processtransitions from one operating point to another.
 13. A method as claimedin claim 8 wherein the step of using a non-linear model includes using arigorous, first principles, non-linear model.
 14. A method as claimed inclaim 8 wherein the subject process is a polymer process.